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Abstract. Among the hundred or so extrasolar planets discovered to date, 19 are orbiting a component of a 
double or multiple star system. In this paper, we discuss the properties of these planets and compare them to the 
characteristics of planets orbiting isolated stars. Although the sample of planets found in multiple star systems is 
not large, some differences between the orbital parameters and the masses of these planets and the ones of planets 
orbiting single st a rs are emerging in the mass-period and in the eccentricity-period diagrams. As pointed out by 
IZucker fc Mazehl J2002f) . the most massive short-period planets are all found in multiple star systems. We show 
here that the planets orbiting in multiple star systems also tend to have a very low eccentricity when their period 
is shorter than about 40 days. These observations seem to indicate that some kind of migration has been at work 
in the history of these systems. The properties of the five short-period planets orbiting in multiple star systems 
seem, however, difficult to explain with the current models of planet formation and evolution, at least if we want 
to invoke a single mechanism to account for all the characteristics of these planets. 
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1. Introduction 

Studies of stellar multiplicity among solar-type stars of the 
solar neighbourhood have shown that about 40% of the G 
and K dwarfs can be cons i dered to be real single stars 
l|Duauennov fc Mavorlll99lUEggenberger et alJl2003|h As 
the majority of solar-type stars belong to double or multi- 
ple star systems 1 , it is of interest to consider the existence 
of planets in such an environment. Searches for extraso- 
lar planets using the radial velocity technique have shown 
that giant planets exist in certain types of multiple star 
systems (see Table 1 for further details). The number of 
such planets is, however, still low, in part because close bi- 
naries are difficult targets for radial velocity surveys and 
were consequently often rejected from the samples. Due to 
the limitations of the available observational techniques, 
most detected objects arc giant (Jupiter-like) planets; the 
existence of smaller mass planets in multiple star systems 
is thus still an open question. 

The orbital characteristics and the mass distribution 
of extrasolar planets can give us an insight into their for- 
mation mechanisms and th eir subsequent evolution. In 
the first paper of this series, lu dry et all (|2003bl : TPaper J) 



discussed the period distribution and the mass-period 
diagram for ex trasolar planets orbitin g single stars. As 
pointed out bv IZucker fe Maze h (2002), planets orbiting 
a component of a multiple star system seem to have dif- 
ferent characteri s tics t han planets orbiting single stars. 
IZucker fc Mazer] l)2002|) showed that there is a significant 
correlation in the mass-period diagram for planets orbit- 
ing single stars, while there may be an anticorrelation in 
this same diagram for planets found in multiple star sys- 
tems. The difference is mainly due to a paucity of massive 
planets with short periods, and to the fact that the most 
massive short-period planets are all found in binaries. 

The characteristics of extrasolar giant planets have 
forced considerable modifications of the standard model of 
planet formation. It is now usually believed that planets 
form within a protoplanetary disc of gas and dust orbiting 
a central star, but the precise modes by which this forma- 
tion tak es place are still debated, especially for giant plan- 
ets (e . g.lPollack et alJl996llBossll9^lBodenheimer et alJ 
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1 In this paper, double and multiple star systems will be 
called multiple star systems. 



l2000l lBossll2000l IWnchterl et alJl2000l: IriosslbOO^ . Two 

major models have been proposed to explain giant planet 
formation (see Sect. 14. each with its advantages and lim- 
itations, but there is currently no consistent model that 
accounts for all the observed characteristics of extrasolar 
planets. Observational constraints are thus needed, not 
only to specify our understanding of planet formation and 
evolution, but also to possibly discriminate between the 
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Table 1. Planets orbiting a component of a multiple star system with confirmed orbital or common proper motion 
(CPM stands for common proper motion and SB for spectroscopic binary). 
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Notes: (a) According to lHalbwachsl l)l98ft) . this pair has only a probability of 60% to be physical. The physical nature of 
this binary has however been confirmed later on the basis of CORAVEL radial velocity measurements (Halbwachs, private 
communication); (b) The multiplicity status of this system has still to be clarified. 
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proposed models. In this context, the detection and the 
characterization of planets orbiting in multiple star sys- 
tems, even if more difficult to carry out than the study of 
planets orbiting isolated stars, may bring new constraints 
and additional information. 

This paper is organized as follows. The sample of plan- 
ets found in multiple star systems is presented in Sect. [5] 
Some trends seen in the statistics are then emphasized in 
Sect.|3J Models of formation and evolution of giant planets 
in binaries are briefly reviewed in Sect. 0] and their pre- 
dictions are compared to the observations in Sect.|S| Our 
conclusions are drawn in Sect . 



2. Known planets in multiple star systems 

Among the extrasolar planets discovered to date, some of 
them arc orbiting a component of a multiple star system. 
Planets have been found around stars known to be part 
of a wide common proper motion pair, known to be in a 
visual binary or in a spectroscopic binary. Alternatively, 
searches for faint companions to stars hosting planets have 
revealed a few new systems. These observations, summa- 
rized in Table 1, show that giant planets can form and 
survive in certain types of multiple star systems. 



3. Statistics of planets in multiple star systems 

Although the sample of planets found in multiple star sys- 
tems is not large, a preliminary comparison between the 
characteristics of these planets and the ones of planets or- 
biting isolated stars can be made. Here, we will discuss the 
mass-period and the eccentricity-period diagrams for ex- 
trasolar planets, focusing on possible differences between 
the two populations. Our sample of planetary candidates 
orbiting a component of a multiple star system consists of 
all the systems listed in Table Q] The total sample of ex- 
trasolar planetary candidates is made of 115 objects 2 with 
a minimum mass A^sini < 18 Mj. The orbital elements 
used for the analysis are the ones deduced from our radial 
velocity data or the most recent version of the ones given 
in the literature. 



3.1. The mass-period diagram 

Figure ^ shows the distribution of all the extrasolar plan- 
etary candidates in the M2sini-logP plane. Two in- 
teresting features emerge from this plot: there are no 
short-period extrasolar planets with a mass M2 sin i > 
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12 3 4 

log(Period) [days] 

Fig. 1. Planetary minimum mass versus orbital period for 
the known extrasolar planetary candidates. Planets orbit- 
ing a single star are represented by open circles, while 
planets orbiting a component of a multiple star system 
are represented by filled circles. The dashed line approx- 
imately delimits the zone where only extrasolar planets 
belonging to multiple star systems are found. 



5Mj, and the most massive short-peri od planets are al- 
most all found in multiple star systems l|Udrv et al . 2002; 
IZucker fc Mazehll2002tfPaper t \. Indeed, planetary candi- 
dates with a mass Mi sin i > 2 M j and a period P < 40 
days are all orbiting a component of a multiple star sys- 
te m, the only e xcept ion being HD 162020 b. As explained 
in lUdrv et ai] l)2002^ . HD 162020b is probably a brown 
dwarf with a true mass much larger than its minimum 
mass, and should therefore be removed from our diagram. 
If it is true that no planet with a mass M2 sin i > 2 Mj 
orbiting a single star has been found with a period P < 40 
days, the two populations of planets are somewhat mixed 
together for periods between 40 and 150 days. The orbital 
period below which there is no more massive planet or- 
biting a single star is thus not well defined, and such a 
unique and well defined limit may, in fact, not exist. 

The paucity of massive short-period planets cannot 
be attributed to observational selection effects since these 
planets are the easiest to detect. Moreover, even if the 
sample of planets orbiting a component of a multiple star 
system is small and incomplete, the presence of a few can- 
didates in a zone of the diagram where there are no other 
planets is significant. We will come back to these differ- 
ences in Sect. 

For periods longer than ^100 days, the distribution 
of the planetary companions in the M2 sin z-log P plane 
is not very different for the two samples. In this period 
range, the mean mass of planets found in multiple star 
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Fig. 2. Eccentricity versus orbital period for all the ex- 
trasolar planetary candidates. Same symbol coding as in 
Fig. ^ The dashed line approximately delimits the zone 
where there are no planets belonging to multiple star sys- 
tems. 



systems is, nevertheless, smaller than the mean mass of 
planets orbiting single stars. This difference comes from 
the fact that no very massive planet has been found on a 
long-period orbit around a component of a multiple star 
system. Again, this cannot be attributed to observational 
selection effects, because several planets with a smaller 
mass and on long-period orbits have been found in multi- 
ple star systems. 

The lack of planets with M 2 sini > 5Mj and P > 100 
days in multiple star systems can, however, be due to the 
small number of planets found in multiple star systems. To 
check this point, we computed the probability that, draw- 
ing a random subsample of 10 planets (i.e. the number of 
long-period planets found in multiple star systems when 
selection effects are taken into account) with P > 100 
days out of the total population of planets with P > 100 
days, we would have no planet with M2sini > 5Mj (see 
Appendix lA.ll for further details). Using the hypergeomet- 
ric distribution ( Appendix lA"| . this probability is 3.25% 
(selection effects have been taken into account by discard- 
ing from the counts the planets with a radial velocity 
semiamplitude K < 15ms _1 ). Thus, if the planets with 
P > 100 days orbiting in multiple star systems have the 
same properties as the planets with P > 100 days orbit- 
ing single stars, the probability not to have a single planet 
with P > 100 days and M2sinz > 5Mj among our sam- 
ple of planets found in multiple star systems is 3.25%. 
This result shows that this trend could be real. It is, how- 
ever, not possible to exclude with a high confidence level 
that the lack of planets with M% sini > 5 Mj and P > 100 
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days may solely be due to small-number statistics and that 
planets with P > 100 days orbiting in multiple star sys- 
tems and around isolated stars may, in fact, belong to the 
same population. A larger sample of long-period planets 
orbiting in multiple star systems will be required to settle 
this question. 

3.2. The eccentricity-period diagram 

The distribution of the extrasolar planetary candidates 
in the e-log P plane is illustrated in Fig. [3 In this dia- 
gram, we note that all the planets with a period P < 40 
days orbiting in multiple star systems have an eccentricity 
smaller than 0.05, whereas longer period planets found in 
multiple star systems can have larger eccentricities. Some 
of the very short-period planets are so close to their par- 
ent star that tidal dissipation in the planet could have 
circ ularized their orbit , even if they were originally eccen- 
tric (|Rasio et al.l ll996'). For longer periods, the orbits are 
not necessarily circularized anymore and any eccentricity 
is possible. For the three planets with a period between 10 
and 40 days orbiting in multiple star systems, the circu- 
larization time (due to tidal dissipation in the planet) is 
r c > 10 12 years. This is clearly too long to explain the low 
eccentricities of these planets by invoking tidal dissipation 
alone. 

The hypergeometric distribution (Appendix [SJ was 
again used to test if the difference observed in this dia- 
gram may solely be due to the small size of the sample 
of short-period planets found in multiple star systems. In 
this case, we computed the probability to have 5 planets 
with logP < 1.6 and e < 0.05 in a subsample of 5 plan- 
ets with logP < 1.6 drawn out of the total population 
of planets with logP < 1.6 (see Appendix IA.2I for more 
details). When selection effects are taken into account (by 
discarding from the counts the planets with a radial ve- 
locity semiamplitude K < 15 ms -1 ), this probability is 
3.77%. This means that the probability to have 5 planets 
with logP < 1.6 and e < 0.05 in the sample of short- 
period planets orbiting in multiple star systems is 3.77% 
if these planets have the same properties as short-period 
planets orbiting isolated stars. On the basis of the current 
samples, there is thus a trend towards a difference in the 
properties of short-period planets orbiting in multiple star 
systems and around isolated stars. The alternative state- 
ment, namely that short-period planets orbiting in multi- 
ple star systems and around isolated stars belong to the 
same population can, nevertheless, not be excluded. The 
observation that the three extrasolar planets with periods 
between 10 and 40 days orbiting a component of a multiple 
star system have very low eccentricities is interesting and 
could be a clue to their formation and/or subsequent evo- 
lution history. We will come back to this point in Sect. El 

For periods longer than ^40 days, there is no signifi- 
cant difference between planets found in multiple star sys- 
tems and around isolated stars. The mean eccentricity is 
similar for the two samples. 



3.3. Remarks 

The limit at ~40 days is valid for both the eccentricity- 
period and the mass-period diagrams, but in the latter 
the limiting period is less well defined. If we plot the evo- 
lution of the mean mass or the highest mass (averaged on 
the three highest values) of planets or biting sin gle stars as 
a function of the period (see Fig. 6 of lPaper lh . there is a 
jump at a period of ^40 days, which reflects the distribu- 
tion of these planets in the mass-period diagram. Thus, 
about the same limiting period is obtained by consider- 
ing the distribution of planets orbiting single stars in the 
mass-period diagram on the one hand, and by consider- 
ing the distribution of planets orbiting in multiple star 
systems in the eccentricity-period diagram on the other 
hand. This is intriguing. 

Even if the orbital parameters of the binaries hosting 
planets are not exactly known, the projected separations 
of these systems (see Table 1) indicate that the five plan- 
ets with a period shorter that 40 days reside in very differ- 
ent types of systems. There is thus no obvious correlation 
between the properties of these planets and the known 
orbital characteristics of the binaries or the star masses. 

The history of planets found in multiple planet systems 
is probably different from the one of "single" planets. For 
the analysis presented here, all planets have been consid- 
ered, but we have checked that our conclusions remain 
unchanged if the planets belonging to multiple planet sys- 
tems are removed from the samples. Among the five short- 
period planets found in multiple star systems, two also 
belong to multiple planet systems. This is something that 
we must keep in mind when discussing the properties of 
these planets. 

4. Models of formation and evolution of giant 
planets in binaries 

Let us now turn to planet formation models and consider 
their predictions regarding the existence and the survival 
of giant planets in binary star systems. We will then be 
in position to compare and confront our observations with 
their results. 

There are different points to take into account when 
considering the formation of giant planets in binaries. 
Indeed, the stellar companion affects all the stages of 
planet formation as well as the subsequent evolution of 
the planet once it has formed. The major points that have 
been studied in the literature are briefly described in this 
section. 

4.1. Binary-disc interactions 

Binary stars can in principle interact with three types of 
discs: two circumstellar and one circumbinary. Transfer of 
angular momentum between the binary and the disc leads 
to a truncation of the inner/outer edge for a circumbi- 
nary/circumstellar disc, respectively. 
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lArtvmowicz fc Lubowl l|l994|) investigated the approx- 
imate sizes of discs as a function of binary mass ratio and 
eccentricity for systems with circumstellar and circumbi- 
nary gaseous discs. Their results show that for a binary 
with a mass parameter /i = M 2 /(M 1 + M 2 ) = 0.3, the 
inner edge of a circumbinary disc is typically located at 
r t — 2.0 <2b (where is the binary semimajor axis) for 
nearly circular binaries and at r t = 3.0 for eb = 0.5. 
The outer edge of a circumprimary disc lies at rt = 0.4 ab 
for nearly circular binaries and at rt = 0.18 ab for eb = 0.5. 
For a circumsecondary disc, the outer truncation radius is 
located near r t = 0.27 for nearly circular binaries and 
near rt = 0.15 ab for eb = 0.4. 

4.2. Long-term stability of orbits 

Assuming that planets can form in binary stars, do long- 
term stabil ity regions exist for pla netary orbits in these 
systems? IHolman fc Wiegertl (^999) studied the long-time 
survival of planets in different regions of phase space near 
a binary star system. Circumprimary as well as circumbi- 
nary planets were studied for different values of the binary 
eccentricity and mass ratio. For a binary with a mass pa- 
rameter (i = 0.3, the largest stable orbit around the pri- 
mary star has a critical semimajor axis r c — 0.37 ab for 
eb = 0.0 or r c = 0.14 ab for eb = 0.5. For the same binary, 
the smallest circumbinary orbit has a critical semimajor 
axis r c = 2.3 Ob for eb = 0.0 or r c = 3.9 ab for eb = 0.5. 

The study of IHolman fc Wiegertl l)l999|) was re- 
stricted to planets i n initi ally circular motion. 
IPilat-Lohinger fc Dvorakl l|2002l) extended this type 
of analysis and determined the variation of the stable 
zone due to an increase of the initial planet eccentricity. 
An increase of the planet eccentricity reduces the stable 
zone, this reduction being of course less pronounced 
that the one due to the same increase of the binary 
eccentricity. 

As shown bv IHolman fc Wiegertl 117) 9!) ) a companion 
star orbiting beyond about 5 times the planetary distance 
is not a serious threat to the long-term stability of plane- 
tary orbits. Nevertheless, this result only applies to orbits 
with a low mutual inclination. Ilnnanen et al 1 J1997I) inves- 
tigated the stability of planetary orbits in binary systems 
with emphasis on the inclination of the orbital planes. As 
an example, they studied the stability of the solar system 
under the presence of an hypothetical distant companion 
placed at 400 AU with differe nt inclinatio ns and masses. 
Due to the Kozai mechanism llKozailll962j) . the system is 
unstable at high inclination when the companion mass is 
larger than 0.05 Mq. A low-mass companion does, how- 
ever, not destabilize the system, even when the inclination 
is high. 

4.3. Giant planet formation in binaries 

Two mechanisms have been proposed to explain giant 
planet formation: core accretion and disc instability (e.g. 



IPollack et al1ll99fit lBosslll997t iMaver et al.ll2002|) . The 

core accretion mechanism begins with the collisional ac- 
cumulation of planetesimals and planetary embryos in a 
protoplanetary disc, as for terrestrial planet formation. In 
the outer parts of the disc, where the amount of solid ma- 
terial is increased by the presence of ices, embryos may 
reach about 10 M e in ~10 6 years and begin to grow an 
atmosphere of disc gas to form giant planets like Jupiter 
and Saturn in ^10 7 years. In the disc instability model, a 
gravitationally unstable disc fragments directly into self- 
gravitating clumps of gas and dust that can contract and 
become giant gaseous protoplanets. Coagulation and sed- 
imentation of dust grains to the protoplanet center could 
form a solid core. This process occurs over a dynamical 
time scale: clump formation and dust grain sedimentation 
proceed nearly simultaneously in ^10 3 years. Though they 
are different, these two mechanisms share a common char- 
acteristic: giant planets should only form in the relatively 
cool outer regions of protoplanetary discs. 

Boss; l)l998l) considered the influence of a binary com- 
panion on giant planet formation via disc instability. 
3D hydrodynamical models of discs with 0.04 Mq were 
evolved in time, subject to the gravity of a binary star 
companion placed on a circular orbit at 40 AU. In the ab- 
sence of the binary companion, the disc is stable, but in 
the presence of the binary companion the disc forms a 
multi- Jupiter-mass protoplanet in 0.002 Myr. 

The evolution of two stars , each o r bited by a circum- 
stellar disc, was simulated bv iNelsonl ((2000) using a two 
dimensional smoothed particle hydrodynamic code. Each 
component of the binary had a mass of 0.5 Mq and a bi- 
nary eccentricity of 0.3 was considered. The system was 
evolved over 2700 yr (8 binary orbits). During and after 
periastron each disc developed strong two-armed spiral 
structures which decayed to a smooth condition over the 
next half binary period; this cycle repeating with little 
variations. The spiral structures decay was due to internal 
heating in the discs, which increased their stability against 
spiral arm growth. Giant planet formation via gravita- 
tional collapse is therefore unlikely in this system. In fact, 
the temperatures in the discs are so high, that some grain 
species, including water ices, are vaporized everywhere. 
Giant planet formation by the core accretion mechanism 
is thus unlikely as well in this system. 

4.4. Evolution of an embedded planet in a binary 

A different approach was considered by iKlev! fcOOlh who 
studied the evolution of a giant planet still embedded in 
a protoplanetary disc around the primary component of 
a binary system. A 1 Mj planet was placed on a circular 
orbit at 5.2 AU from a 1 Mq star. The secondary star had 
a mass of 0.5 Mq and an eccentricity of 0.5. The binary 
semimajor axis was varied from 50 to 100 AU. The sim- 
ulations show that the companion alters the evolutionary 
properties of the planet: the mass accretion rate is in- 
creased and the inward migration time is reduced. 
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4.5. Summary 

The main effect a companion has on a protoplanetary disc 
is a truncation of its radius and an induction of waves 
which, upon dissipation, transfer angular momentum be- 
tween the binary and the disc. Comparing the results pre- 
sented in Sect. l4.ll and in Sect. 14. 2| we see that a planet can 
almost always persists for a long time, wherever it forms 
in a truncated protoplanetary disc. The effect a secondary 
star has on the efficien cy of planet f ormation is, however, 
less clear. According to lNelsor] l)2000j) . the companion has 
a negative influen ce, slowing or inhibiting altogether giant 
planet formation. Boss (1998) claims the opposite, namely 
that giant planet formation via gravitational collapse is 
favoured in binaries. More comprehensive studies will be 
needed, not only to clear up the case of binaries with a sep- 
aration of 40 or 50 AU, but also to explore what happens 
for binaries with different projected separations, mass ra- 
tios and eccentricities. Anyway, we expect the secondary 
star to have an infl uence on pla net formation, at least for 
close binaries, and iKlevI fcOOll) has shown that this re- 
mains true for the subsequent evolution of giant planets 
in binaries. Planets found in multiple stars systems may 
thus have different characteristics than planets orbiting 
isolated stars. 



5. Discussion 

As mentioned in Sect. 14.31 in situ formation is very un- 
likely for short-period Jupiter-mass planets. Formation at 
larger distances followed by inward migration seems to 
be a better explanation to the existence of th ese planets 
(e.g. lLin et alJll99fit iBodenheimer et all EoOOL The high 
masses and low eccentricities of short-period planets or- 
biting in multiple star systems (Sect. [3 and Figs. ^ and 
0) seem also to indicate that some kind of migration has 
been at work in the history of these systems. A few differ- 
ent migration mechanisms have been proposed to explain 
the existence and the characteristics of short-period giant 
planets. We will now briefly discuss two of them and see 
if they might explain some of the features emphasized in 
Sect.H 



5.1. Planet-viscous disc interaction 

One proposed migration mechanism involves the grav- 
itational interaction of a p rotoplanet with the gaseous 
disc out of which it form ed IjGoldreich fc Tremaind ll979. 
Il980t IWardlll986l Il997|) . Subject to such an interaction, 
high mass planets will migrate more slowly than low 
mass planets in a given disc because they create larger 
gaps. Moreover, some of these planets will experience 
mass loss as they come close to the central star. In over- 
all, we thus expect to find more massive planets at in- 
termediate and large semimajor axes, the population of 
close-in objects bein g dominated by smaller mas s planets 
( Trilling et al.ll2002ti . This is indeed observed IjPaner O ) 



and lZucker fc Mazehl l|2002j) have shown that this effect is 
statistically significant. 

Now, if we consider the evolution of a protoplanet or- 
biting the primary star of a binary system, we have seen 
that the presence of the companion alters the evolutionary 
properties of the planet, in partic ular the m igration and 
mass growth rates are enhanced felevll200ll) . These dif- 
ferences may explain why the most massive short-period 
planets are found in multiple star systems: either they are 
more massive because of the higher mass accretion rate, or 
they are massive planets like the ones found with periods 
longer than ^100 days around single stars, but orbiting 
closer-in in multiple star systems because of the higher 
migration rate. Both of these effects are in fact probably 
mixed together and present at the same time. 

Still regarding migration via the gravitational interac- 
tion of a Jupiter-mass protoplanet with a gaseous disc, 
models indicate that a realistic upper limit for the masses 
of closely orbiting giant planets is ^5 M j , if they originate 
in proto planetary discs sim ilar to the minimum-mass solar 
nebula l|Nelson et"aTll2000|) . Examples of large (> 5Mj) 
planets at small orbital distances can, however, be ob- 
tained due to migration in disc s with different masses or 
viscosities l|Trilling et alJll998l) . If such an upper mass 
limit exists for sh o rt-per iod planets, and if the scenario 
proposed by Kiev <)200l[) is correct, we would expect to 
find an upper mass limit for short-period planets orbiting 
in multiple star systems that is larger than the one valid 
for planets orbiting single stars. Therefore, there should 
exist a zone in the mass-period diagram where only plan- 
ets orbiting in multiple star systems are found, the pres- 
ence of a stellar companion being the reason that enables 
a planet to reach a small separation with a mass larger 
than the limit corresponding to planets orbiting isolated 
stars. 



In the simulations bv lKlevlll200l[) . the planet eccentric- 
ity is also modified: it first grows due to the perturbations 
induced by the secondary star, but then declines because 
of the damping action of the disc. The final result is a 
rapid decay of the planet semimajor axis and a damping 
of the initial eccentricity. 

Taken at face value, these arguments may provide an 
explanation for the observation that the most massive 
short-period planets are all found in multiple star systems 
and have very small eccentricities. It should, however, be 
noticed that several of the multiple star systems known 
to host pla nets a r e prob ably very different from the ones 
studied by IKlevI l)200lh . The five planets with a period 
shorter than 40 days orbit in binaries with very differ- 
ent separations (from ~20 to ~1000AU) and it seems not 
likely that the perturbations produced by a wide compan- 
ion would influence the evolution of a protoplanet orbit- 
ing at or below a few AU. This, however, deserves further 
study. 
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5.2. Kozai migration 

Another mechanism that may be at work in wide binaries 
is the so-called Kozai migration Murrav 
In suc h a case, the Kozai mechanism i Kozai 



19621 see also iHolman et all Il997t llnnanen et all Il997t 



Maz eh et al.ll997l) produces large cyclic oscillations of the 



planet eccentricity. During the periods of high eccentric- 
ity, the periastron is small and, consequently, tidal dissi- 
pation becomes important and gradually removes energy 
from the planetary orbit, eventually leading to circulariza- 
tion. The Kozai mechanism coupled with tidal dissipation 
is thus a viable method by which a planet can migrate 
towards the central star. 

Tidal dissipation depends sensitively on the nearest 
approach distance and is important only if the planet can 
reach a high eccentricity. As the eccentricity oscillations 
only depend on the inclination of the planet orbital plane 
relative to the binary orbital plane, the initial inclination 
between these two planes is a key parameter. Moreover, 
for the Kozai mechanism to be effective, the compan- 
ion must provide the dominant contribution to the ap- 
sidal precession of th e planet fsee lHolman et al.1l 19971 and 
IWu fc Murravl EoO^I for more details). The efficiency of 
the Kozai migration is, however, fairly independent of the 
planet mass, and this mechanism will work for planets of 
relatively large masses 

Although the Kozai mechanism may be efficient in bi- 
naries with large semimajor axes, several requirements 
must be simultaneously satisfied for it to operate, and 
such a mechanism will not apply to a large fraction of 
planetary systems. Furthermore, even if the Kozai migra- 
tion has been efficient during a period in the evolution 
of a planet, it does not imply that its orbit is now cir- 
cular. Up to now, the Kozai mechanism and the Kozai 
migration have been considered to explain the high eccen- 
tricity of given planetary candidates such as 16CygBb 
and HP 80606 b llHolman et al.lll997l iMazeh et af1ll997l 
IWu & Murravl2003lK It has never been demonstrated that 
the combination of the Kozai mechanism with tidal dis- 
sipation may account for the existence of close-in planets 
with very low eccentricities. On the other hand, Kozai os- 
cillations are not likely to be currently at work for the 
planets with short semimajor axes, in particular because 
the Kozai mechanism is suppressed by general relativis- 
tic effects. It is thus very unlikely that the low eccen- 
tricity of these planets may be due to the fact that they 
are currently seen in their low-eccentricity phase. It seems 
therefore difficult to explain the characteristics of all the 
short-period planets orbiting in multiple stars systems by 
invoking the Kozai migration alone. 

6. Conclusion 

The characteristics of giant planets found in multiple star 
systems seem to be different from the ones of planets or- 
biting single stars, at least for the short-period planets. 
The major differences are: 



— the most massive (M2 sini > 2Mj) short-period plan- 
ets all orbit in multiple star systems; 

— the planets found in multiple star systems tend to have 
a very low eccentricity when their period is shorter 
than 40 days. 

These observations seem to indicate that migration has 
played an important role in the history of the short-period 
planets orbiting in multiple star systems and that migra- 
tion may be induced differently in binaries than around 
single stars. 

From th e theoretical point of view, it has been shown 
l|Klevll200lh that the presence of a companion star affects 
the properties of a Jupiter-mass planet still embedded in a 
disc around the primary component of a binary by increas- 
ing the migration and mass accretion rates. Alternatively, 
the Kozai mechanism may be (or have been) at work in bi- 
naries hosting planets and will modify some of the orbital 
parameters of the planet. This mechanism can be efficient 
in wide binaries and, coupled with tidal dissipation, it may 
also lead to inward migration. Even if these two mecha- 
nisms may be invoked to explain the characteristics of a 
few planets orbiting in binaries, none of them seem to be 
able to account for all the properties of the five planets 
orbiting in multiple stars systems with a period P < 40 
days. Nonetheless, it is also possible that diverse mecha- 
nisms may have been at work in these systems, but leading 
to a similar final state and similar planet properties. New 
studies dedicated to this issue will be needed to settle this 
question and to find a satisfactory explanation to the ex- 
istence and the characteristics of the short-period planets 
found in multiple star systems. 

From the observational point of view, a larger sample 
of planets orbiting in multiple star systems will be required 
to confirm or refute the preliminary trends emphasized in 
this paper. In this context, the search for planets in mul- 
tiple star systems, even if more difficult to carry out than 
the search for planets around single stars is of importance. 
On the other hand, the characterization of the star sys- 
tems susceptible of hosting planets is underway and could 
bring interesting constraints for the models, thus helping 
our understanding of giant planet formation. 



Appendix A: The hypergeometric distribution 

The hypergeometric distribution models the total number 
of successes in a fixed size sample drawn without replace- 
ment from a finite population of N items of which G are 
labelled success and (N — G) are labelled failure. The hy- 
pergeometric distribution is described by three parame- 
ters: N, the size of the population; G, the total number 
of items with the desired characteristics in the popula- 
tion; and n, the size of the random sample drawn from 
the population. 

The probability distribution of the hypergeometric 
random variable X, the number of successes in a random 
sample of size n selected from the total population is: 
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P(X 



N-G 
n—x 



where x = max(0, n — (N — G)), . . . , min(n, G). This 
probability formula represents the ratio of the number of 
samples containing x successes and (n — x) failures to the 
total number of possible samples of size n. 

To test the statistical significance of the possible dif- 
ference observed in the mass-period or in the eccentricity- 
period diagram, we considered a subsample of planets 
drawn from the total population of extrasolar planets (i.e. 
planets found in multiple star systems and around iso- 
lated stars). A planet was labelled success if it was lo- 
cated within the test zone of the diagram considered and 
failure otherwise. We then computed the probability that 
such a random subsample would give rise to a similar con- 
figuration as the one actually observed for planets found 
in multiple star systems, namely a configuration with the 
same number of planets within the test zone. 

A.l. The mass-period diagram 

We give here more details regarding the statistical signif- 
icance of the lack of planets with M2 sin i > 5 Mj and 
P > 100 days orbiting in multiple star systems (Fig. ^ 
and Sect. 13. 1J) . As explained in Sect. 13.11 the hyperge- 
ometric distribution was used to compute the statisti- 
cal significance of the difference observed. The param- 
eters used were: N — 77, the number of planets with 
P > 100 days and K > 15.0ms" 1 ; G = 21, the num- 
ber of planets with P > 100 days, A^sinz > 5Mj and 
K > 15.0 ms _1 ; n = 10, the number of planets with 
P > 100 days and K > 15.0 ms" 1 found in multiple star 
systems; and x = 0, the number of planets with P > 100 
days, K > 15.0 ms" 1 and M2sini > 5Mj found in mul- 
tiple star systems. Given these parameters, we obtain a 
probability P(X = 0) = 3.25% to have no planet with 
P > 100 days, K > 15.0ms" 1 and M 2 sini > 5 Mj among 
a subsample of 10 planets drawn from the total population 
of planets. 

A. 2. The eccentricity-period diagram 

More details concerning the statistical significance of the 
possible difference observed for short-period planets in 
the eccentricity-period diagram (Fig.0and Sect. 1^.2(1 are 
given here. The parameters used to compute the hyper- 
geometric probability were: N — 25, the number of plan- 
ets with logP < 1.6 and K > 15.0 ms" 1 ; G = 14, the 
number of planets with logP < 1.6, K > 15.0 ms" 1 and 
e < 0.05; n = 5, the number of planets with logP < 1.6 
and K > 15.0 ms" 1 orbiting in multiple star systems; 
and x = 5, the number of planets with logP < 1.6, 
K > 15.0 ms _1 and e < 0.05 found in multiple star sys- 
tems. The probability to have 5 planets with logP < 1.6, 
K > 15.0 ms" 1 and e < 0.05 in a subsample of 5 plan- 



ets drawn from the total population of planets is then 
P(X = 5) = 3.77%. 
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